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Abstract: Air shower universality states that the electromagnetic part of hadron-induced extensive air 
showers (EAS) can be completely described in terms of the primary energy and shower age. In addition, 
simulations show that the muon part is well characterized by an overall normalization which depends on 
the primary particle and hadronic interaction model. We investigate the consequences of EAS universality 
for ground arrays, which sample EAS at large core distances, and show how universality can be used to 
experimentally determine the muon content as well as the primary energy of cosmic ray air showers in a 
model-independent way. 



A ground array detector samples the particles in an 
Extensive Air Shower (EAS) at a limited number 
of points at different distances r from the shower 
axis. From this sample, an observable has to be 
defined to estimate the shower size. To avoid the 
large fluctuations in the signal integrated over all 
distances caused by fluctuations in the shower de- 
velopment, Hillas [1] proposed to use the signal 
at a given distance, S(r), to determine the shower 
size. The optimal distance r op t [2] where experi- 
mental uncertainties in the signal determination are 
minimized is mainly determined by the experiment 
geometry (spacing between ground array detec- 
tors). In this paper, we consider the signal in wa- 
ter Cherenkov detectors as employed by the Auger 
Observatory (r opt = 1000 m). Similar calculations 
can be done for any ground array detectors. 

Using Monte Carlo simulations, S(r op t) is related 
with the energy of the incoming cosmic ray. This 
calibration suffers from large systematics due to 
uncertainties in the hadronic models and the as- 
sumptions that have to be made about the primary 
cosmic ray composition. In this work, we propose 
a new method to determine the calibration in a 
model independent way. Furthermore, this method 
allows us to determine the number of muons pro- 
duced in air showers. 

The method is based on what we will call air 
shower universality [3]: to a remarkable degree of 



precision, EAS can be characterized by only three 
parameters: the primary energy Eq, the depth of 
shower maximum X max , and the overall normal- 
ization of the muon component N^. The param- 
eters X max and Nfj, are linked to the mass of the 
primary particle, ranging from proton to iron, and 
are subject to significant shower-to-shower fluctu- 
ations. All composition and model dependence is 
distilled in these two parameters with clear phys- 
ical interpretation. In addition to determining a 
model-indepedent energy estimator, they can be 
compared with simulations to infer the cosmic ray 
composition and place constraints on hadronic in- 
teraction models. 

Previous studies have demonstrated that the en- 
ergy spectra and angular distributions of electro- 
magnetic particles [4, 5], as well as the lateral dis- 
tribution of energy deposit close to the shower core 
[6] are all universal, i.e. they are functions of 
Eq, X max , and the atmospheric depth X only (the 
dependence on X and X max is commonly put in 
terms of the shower age s). 

By exploiting shower universality, we show that it 
is possible to separate the known shower proper- 
ties, including the electromagnetic particle flux on 
ground and the average depth of shower maximum 
(X max ), from the unknown, the surface detector 
energy scale and the normalization of the muon 
signal at r opt which is tightly correlated with the 
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Figure 1 : Electromagnetic (circles) and muon sig- 
nals (triangles) at r = 1000 m in the shower plane 
vs. distance to ground of the shower maximum, 
for proton (filled symbols) and iron showers (open 
symbols) simulated with QGSJetll/Fluka. 

overall number of muons in a shower. (X max ) as 
a function of energy has been measured with good 
precision by fluorescence detectors, and can also 
be inferred from surface detector variables. 

EAS Universality at large core distances 

In this section we will test shower universality in 
terms of the shower plane signal, i.e. the sig- 
nal generated by particles in a fiducial flat detec- 
tor parallel to the shower plane (orthogonal to the 
shower axis). By avoiding geometric projection ef- 
fects, this allows us to compare showers at differ- 
ent zenith angles. We have assumed a cylindrical 
detector with a top area of 10 m 2 and 1.2 m height 
(similar to the ones used in the Pierre Auger Ob- 
servatory). The response of the detector, simulated 
using Geant 4, is expressed in units of VEM (the 
signal of a vertical, central muon). 

We have generated a library of showers that span 
a zenith angle range of 0° to 70° and an energy 
range of 10 17 eV to 10 20 eV. Showers of pro- 
ton and iron primaries were generated using COR- 
SIKA 6.500/6.502 [7] and the hadronic interaction 
models QGSJetII-03 [8] andFluka [9]. In addition, 
we simulated proton/iron showers at 10 19 eV and 
different zenith angles using other hadronic inter- 
action models (QGSJetII-03/Gheisha2002 [10] and 
Sibyll2.1/Fluka[ll, 12]). 



Figure 2: Electromagnetic signals (in the shower 
plane, r — 1000 m) relative to that of proton- 
QGSJetll. 

The shower-plane signals were separated into sig- 
nals from electromagnetic particles and muons. 
We include the signal from the electromagnetic de- 
cay products of muons (~ 15% of the muon signal) 
in the muon component, the remaining signal be- 
ing the 'pure' electromagnetic component S cm . 

Fig. 1 shows the electromagnetic signal for a core 
distance of 1000 m (circles, proton and iron show- 
ers) as a function of DG = AT groun d — l max , the 
distance from the shower maximum to the detec- 
tor measured along the shower axis (in g/cm 2 ). 
Note that this plot contains showers from all zenith 
angles. Apparently, the signals from proton and 
iron are very similar, though there is a slight shift 
in the overall normalization. This is in violation 
of shower universality, which states that showers 
of the same energy at the same evolutionary stage 
(given by DG) should have the same electromag- 
netic component. 

Fig. 2 shows the electromagnetic signal for dif- 
ferent models and primaries, relative to a refer- 
ence (proton QGSJetll/Fluka). Note that the dif- 
ferent model predictions for a given primary are 
within 5% of each other. There is, however, a sys- 
tematic offset of about 13% between proton and 
iron signals. We also found that the systematic 
differences in the number density of particles are 
smaller, about 8%. This effect persists also when 
comparing signals at the same shower age instead 
of DG. 

Fig. 1 also shows the muon signal (S^, triangles) 
as a function of DG for the same proton and iron 
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Figure 3: Muon ground signals (in the shower 
plane, r = 1000 m) relative to that of proton- 
QGSJetll vs. distance to ground for different pri- 
maries and hadronic interaction models. 
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showers. The dependence of the signal on the pri- 
mary mass (~ 40% between proton and iron) as 
well as the hadronic model is well known. It should 
be stressed that the difference is mostly in the nor- 
malization, not in the functional dependence on 
DG. This is shown clearly by the muon signals 
plotted relative to proton-QGSJetll, Fig. 3. 

Determining the muon normalization 
and energy scale 

The universality of the electromagnetic ground sig- 
nal as well as of the evolution of the muon signal 
can be used to parameterize the total ground sig- 
nal in a model- and primary-independent way. The 
signal at a fixed core distance is then only a func- 
tion of primary energy, distance to ground DG, 
zenith angle 8, and the overall muon normalization. 
The slight primary-dependence of the electromag- 
netic signal enters as a systematic uncertainty in 
the method. Given the measured average depth of 
shower maxmimum (X max ) as a function of en- 
ergy, either from a fluorescence detector (on site 
or a separate experiment) or from ground observ- 
ables, the distance to ground can be directly de- 
termined from the zenith angle for each shower: 
DG = Xq/ cos 8 — (X max ), where Xq stands for 
the vertical depth of the atmosphere at the experi- 
ment site. We parametrize the electromagnetic and 
muon signal as separate Gaisser-Hillas type func- 



Figure 4: Upper panel: the signal parametrization 
Eq. (1) vs. sec 8 for different (black/solid— 1.1, 
red/dashed— 0.5, blue/dotted— 2.0). Lower panel: 
histograms of number of events above the 
parametrized signal in equal exposure bins, ob- 
tained for the same as shown in the upper panel 
from a Monte Carlo data set (see text). 

tions in DG, leaving a normalization factor free for 
the muon signal: 

S(E,8) = S em (E, 8, (X max )) 

+ Nn(E) ■ Sn(0, (X max )) (1) 

Here, denotes a reference muon signal, which 
we take to be proton-QGSJetll at 10 19 eV, and 
N^(E) is the relative muon normalization at this 
energy. Hence, the energy E and N^(E) are the re- 
maining unknowns, which however cannot be dis- 
entangled for individual events in a ground array. 

Fig. 4 (upper panel) shows the zenith angle de- 
pendence of the signal (Eq. (1)) for a fixed en- 
ergy of 10 19 eV and different values of N^. It 
is evident that the smaller the N^, the steeper the 
8 dependence is. We can now use the fact that, 
within statistics, the arrival directions of high en- 
ergy cosmic rays are isotropic. Therefore, we di- 
vide the ground detector data set in equal expo- 
sure bins in zenith angle (sin 2 8 bins). Given a 
muon normalization, we calculate the number of 
events in each bin above a given reference energy 
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(here i? rof =10 9 eV), using Eq. (1). We then adjust 
Nfj,(E rc f) in the signal parametrization Eq. (1) to 
that value which gives an equal number of events 
N(> S(E rc f,8)) in each zenith angle bin (lower 
panel in Fig. 4). For a range of values, we cal- 
culate the x 2 /dof of the event histogram relative to 
a fiat distribution in sin 2 9. This determines the ex- 
perimental value of /V M and its errors. Once is 
determined, Eq. (1) can be used to set the energy 
scale of the experiment. 

In order to prove the feasibility of this method, we 
have simulated 1,000 realizations of a ground ar- 
ray data set with ^2,000 events above 10 19 eV, 
distributed according to the observed cosmic ray 
spectrum and for different primary compositions 
(pure proton, iron, or mixed composition). The 
zenith angle of each shower is sampled from a flat 
sin 2 9 distribution, while X max is obtained from 
the distributions predicted by QGSJetll for each 
primary and energy. is fluctuated according to 
the model predictions. Note that the magnitude of 
fluctuations in X max and are only dependent 
on the primary particle, not the hadronic model. 

Eq. (1) is then used to calculate the signal at 1000 
m from the shower core, 5(1000), which is also 
smeared with an experimental reconstruction accu- 
racy (10% for high signals, and increasing rapidly 
at signals less than 10 VEM). We then applied the 
method described above to calculate the muon nor- 
malization for each simulated data set. We found 
that N/j, is systematically slightly overestimated, 
with the bias mainly depending on composition, 
and only weakly on the detector resolution. For 
pure proton composition, the bias at 10 19 eV was 
found to be around 14% of the true value, while 
for pure iron, it only amounts to a few percent, due 
to the much smaller fluctuations of iron showers. 
This bias can then be subtracted from the deter- 
mined N^. CX p to obtain an estimate of the true N^, 
the uncertainty in the bias entering as an additional 
contribution to the systematic error. Note however 
that a knowledge of (X max ) already places strong 
constraints on the composition. 

Taking into account this knowledge, and assuming 
the observed universality violation and an error of 
(Xmax) of ~ 15g/cm 2 , we found that the total sys- 
tematic uncertainty of achievable is less than 
10%, roughly the statistical error of for this 
data set. 



Conclusions 

Assuming that air shower univerality holds, the 
method presented allows for a measurement of the 
muon content of air showers to better than 10% 
for existing experiments. With similar precision, it 
also determines a converter of signal at ground to 
energy, i.e. a model-independent ground detector 
energy scale. In addition, the measurement of 
can be performed at any energy accessible to the 
experiment. The measured evolution of iV M (E) is a 
further observable of relevance to hadronic models 
and composition. This method has been applied to 
data from the Pierre Auger Observatory [13] yield- 
ing results that constrain hadronic interaction mod- 
els. 
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